Phenomenology of Current-Induced Dynamics in Antiferromagnets 
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We derive a phenomenological theory of current-induced staggered magnetization dynamics in 
antiferromagnets. The theory captures the reactive and dissipative current-induced torques and the 
conventional effects of magnetic fields and damping. A Walker ansatz describes the dc current- 
induced domain-wall motion when there is no dissipation. If magnetic damping and dissipative 
torques are included, the Walker ansatz remains robust when the domain- wall moves slowly. As in 
ferromagnets, the domain-wall velocity is proportional to the ratio between the dissipative-torque 
and the magnetization damping. In addition, a current-driven antiferromagnetic domain-wall ac- 
quires a net magnetic moment. 



In ferromagnets, a spin-polarized current can be used 
to manipulate magnetization via the exchange interac- 
tion. A misalignment between the polarization of the 
current and the local magnetization direction causes a 
spin-transfer torque (STT) on the magnetization because 
of noncollinear spins that precess within the ferromag- 
net. This effect was first theoretically predicted by Slon- 
czewski and Berger [I] and has since garnered abundant 
experimental evidence (for a review, see Ref. The 
STT effect is the reciprocal process of the charge currents 
that are induced by a time- varying magnetic texture [3] . 
A promising commercial application of the STT effect 
utilizes the spin-polarized currents to switch the ferro- 
magnetic layers in spintronic devices, such as in magnetic 
random-access memory or to induce magnetic precession 
for use in high-frequency oscillators in wireless communi- 
cation devices. A limitation of these applications of STT 
is the high levels of critical currents that are required to 
switch the direction of the magnetization. 

Antiferromagnets are ordered spin systems in which 
the magnetic moments of all electrons in each unit cell 
compensate for each other in equilibrium. In spintronic 
devices, antiferromagnets are commonly used to pin fer- 
romagnetic layers via the exchange bias effect. Recent 
theoretical [H [5] and experimental [Bj works indicate, 
however, that current-induced torque effects also appear 
in antiferromagnets. Antiferromagnets also share an- 
other transport property with ferromagnets; namely, the 
anisotropic magnetoresistance (AMR) effect [7] . The an- 
tiferromagnetic AMR effect allows for a detailed experi- 
mental study of the current-induced switching of the an- 
tiferromagnetic layers and the motion of the spatially de- 
pendent antiferromagnetic textures. Therefore, antifer- 
romagnets can replace ferromagnets for use in spintronics 
devices. The applicability of antiferromagnets depends 
on the critical currents that are needed to manipulate 
the staggered magnetization as well as on the resistance 
variation that such a reorientation of the magnetization 
induces. 

Magnetization dynamics in ferromagnets is described 
by the Landau-Lifshitz-Gilbcrt (LLG) equation, which 
has been extended to include STT [8 . In magnetic 



textures, a spin-polarized current contributes two differ- 
ent terms to the LLG equation of motion: a reactive 
torque term and a dissipative torque term The reac- 
tive torque term preserves the macroscopic time-reversal 
invariance of the equation. This term arises from the 
out-of-equilibrium spin density that is induced by drift- 
ing electrons, which have spins that are adiabatically 
guided by the magnetic texture. The dissipative torque 
breaks the time-reversal symmetry and arises from the 
spin-dephasing processes. Except for those systems with 
strong spin-orbit coupling [5] , the magnitude of the dissi- 
pative torque is typically much smaller than the reactive 
torque, but it is not less important [8 . In some fer- 
romagnetic domain-wall systems, the dissipative torque 
dictates the current-induced domain- wall velocity [3]. 

The effect of the reactive STT on the staggered mag- 
netization in antiferromagnets was recently discussed in 
Ref. [5]. The torque has the same physical origin as that 
in ferromagnets. Based on the form of the dissipative 
STT in ferromagnets, Ref. [5] also makes an educated 
guess about the form of the dissipative STT in antiferro- 
magnets. 

In the present paper, we develop a general phenomenol- 
ogy that describes the coupled dynamics of currents and 
the staggered order parameter in isotropic antiferromag- 
nets to the lowest order in spin-texture gradients and pre- 
cession frequency. The antiferromagnet is treated within 
the exchange approximation. For the lowest order in rel- 
ativistic interactions, the exchange forces only depend on 
the relative orientation of the spins. This approximation 
is a good starting point for many conventional ferromag- 
nets and antiferromagnets [10] . including disordered sys- 
tems. In these systems, impurities couple to the spin 
degrees of freedom through random magnetic moments 
or spin-orbit coupling, but impurity averaging restores 
the spin-rotational and sublattice symmetries. We in- 
clude the effects of damping, external magnetic fields, 
and reactive and dissipative torque effects. Our results 
differ from the postulated form of the dissipative torque 
in Ref. [5 , and we explain why. We apply our theory 
to an antiferromagnetic domain- wall system, and find an 
analytic solution in the low current-density regime. Sim- 
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ilar to ferromagnets, we find that the domain-wall veloc- 
ity is proportional to the ratio between the dissipative- 
torque and a bulk damping coefficient. An interesting 
consequence of the current-induced motion is that the 
domain- wall develops a net magnetic moment. Current- 
induced staggered magnetization dynamics can thus be 
observed in two ways: via the AMR effect and via the 
out-of-equilibrium net magnetic moment. 

Our phenomenology is based on the theory of insulat- 
ing antiferromagnets [11] , which is extended to take into 
account the current flow. An important aspect of our 
phenomenology is the exchange approximation that im- 
plies that the total energy is invariant during the simulta- 
neous rotation of all the magnetic moments [llj . Subse- 
quently, when considering the current-induced domain- 
wall motion, we include the magnetic anisotropy phe- 
nomenologically in the free energy, considering that these 
anisotropy energies are very small, e.g., on the scale of 
the critical temperature. 

For clarity, we restrict our treatment to systems in 
which each unit cell in the crystal lattice contains two 
equivalent magnetic sites. In this situation, the antifer- 
romagnet consists of two sublattices with magnetic mo- 
ment densities mi(r,t) and ni2(r, t), such that the total 
magnetization is m(r, t) = mi(r,i) + 1x12(1*, t) and the 
antiferromagnetic order parameter is l(r, f) = mi(r,i) — 
uiafr, t). In equilibrium and in the absence of magnetic 
fields and textures, m vanishes, and 1 is finite and ho- 
mogenous. In the following, we allow the antiferromag- 
net to become distorted into metastable textured states, 
such as domain-walls or vortices, but we require that the 
texture is smooth on the scale of relevant microscopic 
length scales. The texture is parameterized by a slowly 
varying unit vector n(r, t) = l(r, t)/l (I = |l(r, As- 
suming stiff antiferromagnetic ordering, the longitudinal 
dynamics of 1 can be neglected so that the slow dissipa- 
tive dynamics of the system are fully described by the 
directional Neel field n(r, t) along with the transverse 
magnetization m(r,t), which physically corresponds to 
the small relative canting of the magnetic sublattices. 
Constructing the phenomcnological equations of motion 
well below the Neel temperature, we thus impose the 
constraints |n| = 1 and m • n = 0. This starting point is 
analogous to Haldane's mapping of the long-wavelength 
antiferromagnetic action into the nonlinear sigma model 
B2J. 

In addition to rotational invariance, the exchange ap- 
proximation requires that the free energy and the equa- 
tions of motion are invariant under the exchange of the 
two sublattices i.e., that they are invariant under 
the transformations n(r, t) n- — n(r, t) and m(r, t) M- 
m(r,t). The leading-order free energy that satisfies the 



appropriate symmetry requirements is thus 
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Here, we expanded the free energy to the second order in 
the gradients and the magnetization field m(r,i), which 
is coupled to an external magnetic field H. The equations 
of motion for m(r, t) and n(r, t) are found by expanding 
their slow dynamics to the lowest order in the effective 
fields f„ = — 5 n F = An x (V 2 n x n) — m(H • n) and 
f TO = — S m F = —am + n x (H x n). To enforce the 
constraints |n| = 1 and m ■ n = 0, we calculated the 
variational derivatives S m F by varying m normal to a 
fixed n and 5 n F by parallel transporting m on the sphere 
that is parameterized by n. In the absence of electric 
currents, we obtain 113 



h = (jf m - dm) x n, 
rh = (ji n - G 2 n) x n + (ji„ 
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Giih)xm, (3) 



where 7 is the effective gyromagnetic ratio. The dissipa- 
tion power P = rif„+mf m = (Gi/7)m 2 + (G 2 /7)n 2 > 
requires that Gi.,2/7 > 0. The nondissipative equa- 
tions with Gi j2 = are derived in the linearized regime 
in Ref. [TT]. In addition to the appending dissipa- 
tion, we have also added the second term in Eq. (|3j), 
which is quadratic in small deviations from the equi- 
librium, to enforce the constraint m • n = 0. Note 
that such a term naturally appears if one constructs 
the antiferromagnetic equations of motion out of the 
ferromagnetic LLG equations of the constituent mag- 
netic sublattices. Eq. ([2])-([3]) can be reduced to a sin- 
gle equation for the Neel field (without dissipation): 

n x ri = 7 2 an x AV 2 n — H (H ■ n) /a + H x 11/70 — 

27 (H ■ n) n. This equation agrees with the equation 
that is derived in Ref. [TUj from the Lagrangian density 
C = (ri/7 - H x n) 2 /2a - A (Vn) 2 /2. 

We also make use of the linearized equations in the 
Landau-Lifshitz form: 

n = 7(f m x n + Gif„) , m = 7(f„ x n + G 2 f m ) , (4) 

where 7 = 7/ (1 + G1G2). The Onsager reciprocity re- 
lations between the two fields require the gyromagnetic 
ratios to be the same in the two equations (see below). 

The equations ^ and Q describe the evolution of 
an electrically open antiferromagnet. Next, we include 
the effect of itinerant electrons on the long-wavelength 
dynamics of n(r, t) and m(r, t) by adding torque terms 
that arise from the currents that are induced by an ex- 
ternal electric field E. To this end, we are guided by 
the rotational symmetry requirements and the Onsager 
reciprocity relations. The Onsager reciprocity relations 
apply to a system that is described by several param- 
eters {qi\i = 1, . . . , N} for which the rate of change <ji 
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is induced by the thermodynamic forces fa = —d qi F, 
and state that the off-diagonal linear response coeffi- 
cients in the equations = J2f=i^ijfj are related by 
Ly(H,M) = e^L^-H.-M), where e, = 1 ( Ci = -1) 
if qi is even (odd) under time reversal. Here, M repre- 
sents any possible equilibrium magnetic order. The fields 
that describe the collective magnetic dynamics in the an- 
tiferromagnet are n(r,t) and m(r,i), and the associated 
conjugate forces are f„ and f m , respectively. In the diffu- 
sive regime, the charge transferred by the current-density 
j is conjugate to the electric field such that i q = E. The 
response coefficients that are needed are the response 
matrices L n ,q and L m . q , which describe the dynamics 
of n and m that are induced by the electric field. Be- 
cause the magnetization is odd and the charge is even 
under the time reversal, Onsager's theorem implies that 
ini(m»),9j( n > m ) = -Ajj.nitm^C-n, -m), where L q ,„(m) 
describes the charge currents that are pumped by f„ (f m ), 
and the (±n, ±m) arguments denote an equilibrium tex- 
ture. 

To derive the STT terms, it is convenient to be- 
gin by phenomenologically constructing the magnetically 
pumped charge current jP um P, which yields £ q , n (m)i an d 
then invoke Onsager's theorem to obtain L n ( m ) jq . For the 
lowest order of the space-time gradients and the magne- 
tization field m, we can write three pumping terms that 
satisfy the appropriate exchange and spatial symmetries: 
n(m x dia), n-c^n, andn-(n x dim). However, because 
the last term is quadratic in the small deviations from an 
equilibrium state (in the absence of magnetic fields), we 
disregard it in the following. Thus, the leading-order phe- 
nomenological expression for the pumped charge current 
is as follows: 
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pump / 



j a = rjn ■ (m x c^n) + j3 n • diii 



= 7 [fa + GiP)din ■ f„ + [fi - G 2 v)n x 9,11 • f m ] , (5) 

where we have utilized Eq. Q and have scaled the cur- 
rent density with the conductivity a. Here, rj ((3) is a 
phenomenological parameter. Later, it becomes clear 
that r\ (/3) parameterizes the adiabatic (non-adiabatic) 
torque because the term is even (odd) under time re- 
versal. Eq. (|5| yields the response coefficients L qun . = 
aj(?]+Gif3)d i n j and L qu , n . = a^{t3-G 2 -n){nxd i n) r Us- 
ing the Onsager reciprocity relations and Ohm's law for 
the drift current (j = ctE), leads to the STT terms t„ = 
j(t]+GiP)(3- V)n and r m = -7(/3-G 2 ?y)nx (j- V)n for 
the Neel and the magnetization field, respectively, which 
are added on the right side of the equations of motion 
in Eq. Q. Transforming these torques back to the LLG 
form of the equations, i.e., Eqs. Q and ([3]), yields the 
following: 

n = (jf m - Girh) x n + i]-f(j ■ V)n , (6) 
m = [ 7 f„ - G 2 h + /3 7 {j • V) n] x n + r nl . (7) 

r n i = (-f{ m - Girii) x m - 777 [m • ( j ■ V)n] n are the 
simplest nonlinear terms that arc added here to enforce 



the constraint m • n = 0. We disregard such higher-order 
terms in the following. 

From now on we use the simplified notation for the 
effective forces [T3], which allows to more readily solve 
for m in terms of n. Combining Eqs. ^ and Q, we see 
that the magnetization field is fully determined by the 
order parameter n and its dynamics: 
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n x H + -h — Gttn - (77 + GM ■ V)n 

7 



(8) 

Substituting this into Eq. Q allows us to derive a closed 
equation for the Neel field to the linear order in the out- 
of-equilibrium deviations m, <9 t n, j, and H [13] : 

fi/7 = - n x H + Gif„ + (77 + GtflO ■ V)n 

+ a [ 7 f„ - G 2 h + 1 p{ 3 ■ V)n] . (9) 

Eq. ([5])-(|9]) are our main results, which describe a gen- 
eral phenomenological theory of weakly excited current- 
induced dynamics in conducting antiferromagnets and 
charge pumping that arises from moving textures. The 
reactive torque in Eq. ([9]), which is proportional to 
(77 + Gi/3), was first found in Ref. [5]. The dissipative 
torque term, which is proportional to a/3, and the ef- 
fects of magnetization damping are new terms that have 
not been derived before. The consideration of the charge 
pumping that occurs when moving antifcrromagnetic tex- 
tures is also new. Ref. [5] suggests a dissipative torque of 
the form /3n x (j • V)n in Eq. Q. This term breaks the 
n i— > — n invariance of the equation and therefore can- 
not appear in the exchange approximation of equivalent 
magnetic sublattices. 

As an application of our theory, we consider an anti- 
ferromagnetic domain- wall system and study the current- 
induced domain-wall motion. For clarity, we set the ex- 
ternal magnetic field to equal zero from this point of the 
paper on. Domain-walls can be created in systems with 
anisotropy, which is added phenomenologically to the free 
energy as follows: F [m, n] — > F[m, n] + W^[n], where 
W [n] = J dr (Kj_riy/2 — K z n 2 z /2) is the anisotropy en- 
ergy (K± , K z > 0) . A local minima of the above en- 
ergy functional is a Neel wall that rotates in the xz 
plane, where the local magnetization direction is n x = 
1/ cosh(z/A tt ,) and n z = tanh(z/A w ). A^, = A/K z is 
the domain- wall width. This equilibrium domain- wall 
texture is denoted by no(r,t) below. 

In the following, we study how the domain-wall moves 
in response to a current along the z axis. Let us first con- 
sider the case of no magnetization damping and (3 = 0. 
In this case, it follows from Eq. (|6])-([7]) that m(r, t) = 
and n(r, t) = no(z — r w (t)) is an exact solution of the 
equations, with the domain-wall velocity r w = —jrjj. 
By including the magnetization damping and the dis- 
sipative torque, a local magnetic moment density devel- 
ops. The torque T m = —jj ((3 ~ G 2 ?7) n x d z n induces 
a magnetic moment density along the y-axis that should 
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eventually approach a finite value due to the opposite 
acting damping term G2f m in Eq. Q. Thus, to find 
a stationary solution for Eq. ([6|-([7|, we use the ansatz 
n(r,t) — n (z—r w (t)) and m(r,t) — m (i)nx9 z ii. Here, 
mo (t) parameterizes the magnitude of the local magnetic 
moment density. Substituting these two expressions into 
Eq. (J©])- ([t]) , produces the following equations for the two 
parameters r w {t) and mo(t): 



f w = 7 [am - (77 + Gi/3) j] , 
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(G 2 r)-P)j-G 2 a 1 + ^ 



jG 2 aX v 



(10) 
m . (11) 



(ID 



The second term inside the last parenthesis in Eq 
is position-dependent because of n z . When this term 
is negligible, the ansatz becomes a good approximation 
of Eq. ( 10 )- ( 1 1 1 . This low current-density regime cor- 
responds to systems for which the characteristic intrin- 
sic relaxation time (^yG^a) -1 of the antifcrromagnetic 
system is much smaller than the timescale X w /r w ; i.e., 
the domain-wall moves a small distance as compared 
to the domain-wall width during the relaxation time. 
In this regime, m approaches a finite stable value of 
mo = — (0 — G 2 rf) j/ (G 2 a), and the domain-wall moves 
at a constant velocity r w = —^f3]/G 2 . Similarly to that 
for ferromagnets, the velocity is proportional to the ra- 
tio between the dissipative torque and a Gilbert damp- 
ing coefficient. In contrast to ferromagnets, this result 
is independent of the uniaxial anisotropy K± such that 
the Walker ansatz has a much wider range of applicabil- 
ity. An interesting consequence of the current-induced 
domain-wall motion is that the system develops a finite 
magnetic moment in the domain-wall region. This effect 
may be an alternative to the AMR effect for the mea- 
surement of domain- wall motion. 

To verify that the system approaches the above sta- 
tionary solution in the relevant regime, we conducted 
a micromagnetic simulation of a one-dimensional sys- 
tem based on Eq. ([6])-([7|. For the numerical calcula- 
tion, we wrote the equations in a dimcnsionless form 
by scaling the z axis with the lattice constant ai and 
the time axis with (^A*) -1 . Here, A* = A/ (laf). 
We considered a domain-wall system with a domain- 
wall width of Au, = 20a/. The anisotropy and damp- 
ing parameters are all A* = 10, K Z /(IA*) = (20)~ 2 , 
K±/(IA*) = 0.1 and Gxl = G 2 /l = 0.01. For the STT 
torque parameters we used j (77 + Gi/J) / (a/A*) = 0.1 
and j(J3- G 2 ri) / (a t A*l) = 0.001. With these values, 
the stationary solution implies that = mo/ (ail) will 
approach 0.01 and that the wall will move a distance of 
2ai during the relaxation time (^yG^a) -1 . We therefore 
expect the system to be in the relevant regime for which 
the stationary solution is valid. 

Fig.[T]shows the micromagnetic simulation of the above 
system. We see that the domain-wall velocity follows 
Eq. ( |10[ ) nearly perfectly and that velocity and mo ap- 
proach the expected stationary values. It should be noted 
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FIG. 1: The blue line shows the domain- wall velocity, which 
was determined using a micromagnetic simulation, as a func- 
tion of time when a current is applied at t = 0. The ve- 
locity follows the analytic expression in Eq 
proaches the stationary value r w = —"ffij/Gi. Inset 
time evolution of mo. The parameter approaches the sta- 
tionary value mo = — (j3 — G2f)) 3/ (G2«). All of the re- 
sults are given in dimensionless quantities, a* = al/A*, and 
P* =j(v + G 1 /3)/(a l A*). 



that our ansatz breaks down when (7G2&) -1 is not much 
smaller than \ w /r w ; however, a more detailed study of 
this regime is beyond the scope of this manuscript. 

The typical values of the /3/G-ratio in antiferromagnets 
is an interesting issue for future experiments. This ratio 
can be probed by measuring the domain-wall velocity as 
a function of the current density or by measuring the 
reciprocal process, which is voltage that is induced by a 
moving domain- wall. In the last case, one could initiate a 
domain-wall motion using a short current pulse and then 
measure the voltage echo from the reciprocal process. 

In conclusion, we have derived a general phenomeno- 
logical theory of current-induced dynamics in antiferro- 
magnets and have applied the theory to the study of 
current-induced domain-wall motion. We found that the 
domain-wall developed a net magnetic moment during 
the current-induced motion and that the domain- wall ve- 
locity was proportional to the ratio between the dissipa- 
tive torque parameter and a damping parameter. 
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